Ab-initio prediction of magnetoelectricity 
in infinite-layer iron oxides 
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Density functional based simulations are employed to explore magnetoelectric effects in iron-based 
oxides, showing a unique layered structure. We theoretically predict CaFe02 to be a promising 
magnetoelectric, showing magnetically-controlled large electric polarization, possibly even above 
room temperature. The cross coupling between magnetic and dipolar degrees of freedom needs, 
as main ingredients, Fe-site spin-orbit coupling and a spin-dependent O p - Fe d hybridization, 
along with structural constraints related to the non-centrosymmetric point group and the peculiar 
geometry characterized by "flattened" Fe04 tetrahedrons. In order to enhance magnetoelectric 
effects, we performed a materials-design leading to a novel and optimized system, MgFe02, where 
the larger O4 tetrahedral distortion leads to a stronger polarization. 

PACS numbers: Valid PACS appear here 
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Introduction. — Aiming at the mutual control of elec- 
tricity and magnetism, the magnetoelectric (ME) cou- 
pling between the magnetic order parameter (A/) and the 
ferroelectric one (P) has been well studied in recent years 
in the context of multiferroic oxides. [l| Microscopically, 
common mechanisms for ME effects can be classified into 
either the inverse Dzyaloshinskii-Moriya mechanism [2| 
(or the equivalent spin-current mechanism [3|) showing 
P ex J2ij ^ij ^ i'^i ^ •^j) between neighboring spins 
connected by a vector e^, or the inverse Goodcnough- 
Kanamori (or exchange-striction) mechanism |J| show- 
ing P oc J2ij ■^iji'^i ' '^j) with exchange integral Jy. 
Recently, an alternative mechanism, denoted as "spin- 
dependent p-d hybridization" was proposed. |5|,|6| There, 
the spin-orbit-coupling (SOC) "asymmetrizes" the p-d 
hybridization between the transition metal (TM) and 
the surrounding ligands, inducing an electric polariza- 
tion P ex y^-(S'; • e')^e'-, where e'- labels the vectors 
connecting the TM to the ligand ions, which explains 
ferroclcctricity in CuFc02. It has been more recently 
reported that the mechanism can be responsible for the 
polarization observed in Ba2CoGe207 (BCGO), where 
two neighboring Co spins are aligned in an antifcrromag- 
netic (AFM) configuration, [7| as confirmed by our previ- 
ous DFT study. [8[ In fact, in the C0O4 tetrahedral coor- 
dination, the occupied non-bonding x^ — y^ state for mi- 
nority spins is mixed with unoccupied bonding yz and zx 
states through SOC, leading to asymmetric pd hybridiza- 
tion and, in turn, to a net polarization. Incidentally, we 
found a relevant ME mechanism in magnetite, Fe304, 
where the complicated charge-ordering pattern with po- 
lar Cc space group causes a large spontaneous polariza- 
tion, whereas SOC causes a small spin-dependent change 
in the polarization under the ferrimagnetic order. [9| 

Guided by the same theoretical analysis performed for 
BCGO and Fe3 04, in the search for novel ME materials. 



we identified the requirements for the spin-dependent pd 
hybridization mechanism as follows: i) a non-polar but 
non-centrosymmetric point symmetry at magnetic sites, 
whose symmetry is lowered into polar group under a spec- 
ified magnetic order; ii) weak magnetic anisotropy (so 
that spins easily follow the magnetic field), at least in a 
plane, iii) strong pd hybridization and strong SOC needed 
to achieve a large ME effect. Whereas the third condition 
can be fulfilled by using either 4d or bd elements instead 
of 3d transition metals, our approach consists in explor- 
ing different crystal structures to enhance the asymmetric 
pd hybridization to realize larger polarization. 

Our strategy is rather simple. We searched for ox- 
ides which have space group similar to BCGO and iden- 
tified CaFc02 as a promising candidate. Indeed, CaFe02 
crystalizcs in non-centrosymmetric and non-polar PA2im 
(#113) tetragonal layered structure. The crystal struc- 
ture shows peculiar Fe04 squares heavily distorted to 
tetrahedrons (Fig. [1] (a)). Experimentally it can be 
synthesized by removing only apical oxygens from FeOg 
structure in precursor Brownmillerite CaFe02.5 by using 
reductive hydride CaH2. The peculiar Fe04 tetrahedron 
is unique to CaFc02. For example, no distortion has been 
reported in the infinite-layer cuprate ACUO2 (A= Ca, 
Sr, Ba), well studied in the context of high-temperature 
superconductivity [10|, as well as SrFe02, having high- 
spin Fe ions in planar square coordination. [ll| So far, the 
recently synthesized CaFe02 uniquely shows the corruga- 
tion of Fe02 planes, where the O ion is slightly deviated 
by O.32A from the ab plane. |12| 

We note that CaFe02 exhibits G-type antifcrromag- 
netic order (both in-plane and inter-plane AFM coupling) 
below a remarkably high Neel temperature, Tn=420K — 
even above room temperature — . Four Fe atoms in the 
unit ceh, Fel at (0, 0, 0), Fe2 at (1/2, 1/2, 0), Fe3 at 
(0, 0, 1/2), and Fe4 at (1/2, 1/2, 1/2), are responsible 
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FIG. 1: (a) Crystal structure and (b) projected image in the 
ab plane of CaFe02 with the P-42i7Ti space group; Fe ions 
(located in flat O4 tetrahedrons) lie in c=0 planes whereas 
Ca ion lies in c—1/2 planes, (c) Fe spin configurations: (i), 
(ii) coUinear AFM and (iii) non-coUinear spin-canted under 
applied ii"//l-10. Note that the AFM inter-layer coupling is 
not shown here, (d) Development of MO4 tetrahedron with 
respect to the planar square coordination by substituting M 
by Sr, Ca and Mg. 



for the magnetism (the magnetic unit cell is double along 
the c axis compared to the structural unit cell). 

In this study, by means of DFT calculations, we first in- 
vestigate the steric effect of A site substitution (A= Mg, 
Ca, Sr, Ba) in AFe02 and the stability of the infinite- 
layer structure; then, we investigate the ME effect in- 
duced by varying the magnetic order under SOC. 

Symmetry Analysis. — Although the space group of 
crystal structure of CaFe02 is the same as BCGO, the G- 
type AFM order may give a different ME behavior com- 
pared to the C-type (only in-plane AFM coupling) AFM 
order in BCGO. In order to characterize possible ME ef- 
fects, we first briefly perform a group theory analysis and 
discuss its implications in the framework of Landau the- 
ory of phase transitions. [13J In the parent P42iml' space 
group with eight symmetry operations {E, C2(z), 2S'4, 
2C'2(a:,y), 2ad\ plus time-reversal {1'}, the magnetic order 
leads to a lowered symmetry. We define the order param- 
eters F = Si + S2 + S3 + Si'ATLAA^ S1-S2-SS. + Si 
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FIG. 2: (a) DFT results for Pc as a function of the collinear 
spin angle (j) in the ab plane, fitted to sine curves. Spin con- 
figurations in the ab plane shown by arrows, (b) DFT results 
for Pc as a function of the noncollinear spin-canting angle (f)' 
in the ab plane, fitted to cosine curves (solid line). 
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as the ferromagnetic (FM) and G-type AFM combina- 
tion of Fel..4 spins, respectively. Since the G-type AFM 
order doubles the unit cell along the c direction, we have 
to take into account the c translation {c =-|-(001)} in 
this analysis, (i.e. AFM order is not invariant under the 
c translation.) 

Using the transformation rules given in Table HI the 
thermodynamic free energies in terms of all the possible 
ME coupling terms (of the form P ■ Af ^ . which arc in- 
variant under any symmetry operations) results in the 
simple exression: 



FmE = CAPcAaAb + CpPcFaFb 



(1) 



The dielectric energy is written as Fde = —P'^I'^X- Here, 
caiCf and x (henceforth set as 1) are constants. P is 
then evaluated at the minima of i^ = -Fme + ^Fde, reading 



Pa^Ph= 0, 



Pc^CA^aAb+Cpi^a-Fb 



(2) 



We observe that this result means that P^ can origi- 
nate from either FM or AFM order, but not from their 



TABLE I; Matrices of the generators of space group P42iml' 
in the representations spanned by F, A, and P. The group 
elements denote the identity, 7r-rotation, 7r/2-rotoinversion 
5^(=7C^), screw C2i/ + (||0), time-reversal 1' and transla- 
tion c. Labels of irreducible representation (IR) are taken 
from the ISODISTORT program, [li] 
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combination. 

From now on, wc will focus on the Pc behavior assum- 
ing a canted AFM configuration under an applied mag- 
netic field, i.e., we first simultaneously counter-clock- wise 
rotate four Fe spins in the ah plane keeping the AFM 
order with the angle 4> from the a axis, then we cant 
spins by an angle </>', as depicted in Fig. [TJ Accordingly, 
we set 51 = 54 = S'(cos(0 -I- 0'), sin(0 -|- 0'), 0) and 
82 = Sz ^ S{ — cos((/) — (/)'), — sin(0 — 0'), 0), ending up 
with 



Pc = 405*2 sin 20 (cos 20' -t- h) , 



(3) 



where a = ca + cp and 6 = ca + cp — 2. 

A spontaneous Pc is therefore induced in the Ano 
or Ai_io ((/)=±45°) order but not in the Aioo (Aqio) 
one. Analogously to previous studies [8|, l9(, the non- 
magnetic group lacks inversion symmetry, but the sym- 
metries which prohibit Pc (e.g. C2y rotation) are broken 
by the Ano magnetic order. Finally, Eq. ([3]) at fixed 
(0') gives the simple 0'(0)-dependence, Pdfj)) oc sin 20 
and Pcifj)') ^ cos 20'. This result shows a slightly differ- 
ent aspect from the P induced in C-AFM BCGO0,[3, 
where Pf ^°°(0') oc cos(20'-/?) and the phase shift /3 de- 
pends on the non-zero cap coefficient. Notably, the phase 
shift of the cosine curve, /3, is lost in G-AFM CaFe02. 
Physically, a close inspection of the magnetic and ionic 
configuration shows that the cancellation is due to the C 
type interlayer AFM coupling and the consequent spin 
canting pattern in different layers (see Fig. [Ijc)-(iii)). 

DFT analysis — In order to quantitatively confirm 
the ME behavior, we performed DFT calculations on 
AFe02 {A = Mg, Ca, Sr, Ba) by using the VASP code[l3 
with GGA-PBE-HC/ potential (C/=4 eV for Fe-d state, 
taken from Ref.[12J). Note that bare GGA calcula- 
tions result in quasi-metallic state with small gap, e.g. 
Eg <0.1eVin MgFe02. After structural optimization un- 
der G-AFM configuration without SOC, we found that 
MgFe02 and CaFe02 show a corrugation in the Fe02 
layer in non-centrosymmetric P-42ito structure, whereas 



TABLE II: Magnetic anisotropy energy (MAE) (meV/Fe) 
obtained by comparing the total energy with different spin 
directions under SOC and for different values of U in the 
GGA-|-[/ scheme. Spin and orbital moment (S and L) (/^b) 
are also reported for S / /{IQQ). 
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SrFe02 and BaFe02 show a "ffat" layer in centrosym- 
metric PA/mmm structure. In the Fe04 tetrahedra of 
MgFe02 and CaFe02, the JT-active Fe'^+ (dP) ion shows 



e1^t2ge^^t2g occupied orbital states, where the almost- 
non-bonding Sz^ — r" 



r^ states perpendicular to the ab plane 
is the lowest energy orbital state. (Here we choose a local 
frame x//a, y//b, z//c and neglect the small tilting of 
O4 tetrahedra with respect to the crystalline axis.) The 
calculated small magnetic anisotropy in the easy ab plane 
and hard c axis (Tbl. HIl grant an easy control over the 
spins by applied magnetic field, H; even under a small 
magnetic field, the spins flop to be perpendicular to H 
and then cant in order to reduce the Zeeman energy. 

Imposing the coUinear AFM configuration, we simul- 
taneously rotate the Fe spins in the ab plane. We eval- 
uated the ME effect as the change of P (calculated by 
Berry phase method [16[) induced by the rotation of M 
with respect to the crystalline axes, including the opti- 
mization of the atomic coordinates. Figure [2] (a) shows 
Pc as a function of the spin-rotation angle 0: it is consis- 
tent with the previously discussed Landau analysis, being 
Pc oc sin 20. As summarized in Tbl lllH the purely elec- 
tronic contribution via SOC at fixed atomic structure, 
Pcicc=13/iC/m2, is strongly enhanced (up to 62fiG/in^) 
when internal atomic coordinates are optimized. This sit- 
uation is similar to previously reported calculations on 
multiferroic TbMnOs, where the purely electronic con- 
tribution P'=''='^=32/iC/m2 is enhanced up to peiec-hion^_ 
4:67 fj,C/m^ by ionic relaxation. [17f In CaFe02, the max- 
imum value of the calculated polarization is 62fiC/m^, 
slightly larger than the corresponding value in BCGO 
(P=47/xC/m2). Furthermore, MgFe02 is predicted to 
show a much larger polarization, P=327fiC/m^, compa- 
rable with TbMnOa. 

TABLE III: Deviation of Fe-O-Fe bond angle from square 
configuration and Pc under Ano magnetic ordering, calcu- 
lated in the fixed non-polar atomic structure (-Poxp) and with 
optimized polar structure (Popt). Experimentally, the Pc in 
Alio Ba2CoGe207 is about 100/iC/m^[3j, whereas polariza- 
tion in CaFe02 and MgFe02 has not been measured. 
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The strong increase of P in MgFe02 is attributed to 
the larger O4 tetrahedral distortion with respect to the 
centrosymmetric flat layer structure. In fact, we also 
performed calculations on BeFe02, showing even more 
distorted tetrahedrons. Whereas we estimated a signif- 
icantly large polarization, P ^2000 nC/rri^, the chemical 
stability of the crystal structure with such a small ionic 
radius ion is questionable. It is shown in Tbl lllll that 
the ionic contribution (mostly oxygen ionic displacement) 
plays a crucial role in total P. The structural difference 
of iron-based oxides compared to BCGO is that Fe04 
tetrahedra are corner-shared, whereas C0O4 tetrahedra 
are intercalated by Ge04 tetrahedra in BCGO. This leads 
to a cooperative effect between neighboring Fe04 tetra- 
hedra, so as to enhance the ionic displacement of O ion 
shared by two Fe04 tetrahedra. 

The ionic displacement driven under Ano magnetic or- 
der in CaFe02 is characterized by two phonon modes, 
Fi and T3. While the Fi phonon mode doesn't change 
the original space group, the F3 mode reduces the sparse 
group from P-42iml' to polar Cmm21' . Under the F3 
phonon mode, MgFe02 exhibits Mg and O ionic displace- 
ments: both ions are shifted along the z direction, by 
0.0006 A and by 0.0011 A, respectively (displacements 
characterized by Al(l) and ^2(0)). This means that 
ferroelectricity is strongly coupled, through magnetism, 
with lattice distortions in a sort of magneto strictive- 
piezoelectric effect. 

We look then at the spin-canting effect induced by an 
applied field iliio- Figure [2] (b) shows the change in 
Pc induced by artificially canting the spins by an an- 
gle (/)', starting from the Ai_io AFM configuration. In 
agreement with the Landau theory analysis, Pc evolves 
as cos 2(0') -I- const.. 

Single- site SOC induced ME effect — The microscopic 
origin of P in CaFe02 can be explained in terms of spin- 
orbit coupling term treated as second order perturbation. 
In the flat tetrahedral Fe04 coordination (Fe is at 2a site 
with —4.. symmetry), Fe-d electrons occupy 3z^ — r^, 
x^ — y^, degenerate yz/ zx, and xy orbital states in order, 
which is consistent with order of the nonbonding {i.e. 
spatially avoiding O ligands) character. [19j Comparing 
the DOS of CaFe02 and MgFe02, the electronic struc- 
tures are very similar, although the crystal field splitting 
between 82;^ — r^ and x^ — y^ (=Ai) is comparatively 
weakened in more tetragonally distorted MgFe02. 

The spin orbit coupling (SOC) term is described as 
Hsoc = A J2a,a'i^\^ • S\a')d'ld^,, where the matrix 
elements can be expressed as a function of the polar 
and azimuthal angles {9, (p) defining a local reference 
for the spin-quantization axis[l8|. SOC-induced mixing 
of the local d levels between occupied and unoccupied 
states lifts the degeneracies of yz and zx manifolds 
and implies asymmetric hybridizations with the ligand 
oxygens, that ultimately induces a local dipole moment 
along the c axis. Large SOC mixing is expected at Ai 
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FIG. 3: Density of states (DOS) around the Fermi energy 
(£=0) of CaFe02 and MgFe02. The d-orbital projected DOS 
and SOC mixing between occupied and unoccupied states are 
labeled. 
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FIG. 4: (a) Schematic picture of unoccupied Fe-d orbital state 
which is mixed into occupied state under SOC. The orbital 
shape changes by rotating Fe spin, (b) Charge density iso- 
surface of the highest occupied band (i.e. Fe-d iz^-r'^ orbital 
state) in MgFe02 under Ano magnetic order. Spin direction 
is shown by an orange arrow, (c) Change in the charge den- 
sity [5p) by rotating Fe spin from [110] to [1-10] direction 
(positive shown in blue color; negative, in red). Direction of 
polarization caused by the charge difference (±5) is shown by 
an open arrow. SOC is enhanced by factor of 10. 



between [Sz^ — r^|) and |a;^— y^J,), A2 between \xy^) 
and [a:^ — y^4) and A3 between \xy\) and jyZjZxJ,), 
shown in FiglS] Among them, Ai and A3 are relevant 
to ferroelectricity. Under this SOC-related mixing, the 
occupied states are modified by small contribution from 
the unoccupied states; 



fects, with a giant polarization tuned by magnetic fields. 



S\3z^-r^\) =-^{cos<j)\yzi)-sm(l)\zxi)) (4) 
Slxyt) = -2^icos(l}\yzi)+sm(t)\zxi)), (5) 

where is already set as 90° . Comparing the coefficients 
in above equations |3][n] and considering Ai and A3 to 
have a similar magnitude (see FiglS]), we observe that 
the SOC mixing involving Ai contributes to the change 
in orbital occupancy from the change in spin azimuth an- 
gle 4> by an amount which is about 3 times larger than the 
SOC mixing involving A3. Therefore, the energetically 
highest occupied state |3z^ — r^|) is mixed with unoccu- 
pied (cos (j) \yzi) + sm(j)\zxl)) orbital states, which mod- 
ifies the shape according to the spin rotation (see Fig|3] 
(a)). This gives rise to the asymmetric pd hybridiza- 
tion, enhancing the bonding character with upper oxy- 
gen states (when S//[110]) or with lower oxygen states 
(when S//[l-10]). The change in the charge density of 
the occupied state by spin rotation indeed exhibits an 
asymmetric pd hybridization (shown in Fig. |4] (c)). In 
this way, a local dipole, Pc oc sin2(/), develops along c, as 
predicted from the functional form, P oc ^ ■ - {Si ■ e^)^ e'j. 
[5|-l7| Such a mixing of local d-levels nicely explains why 
the Pc size is maximum at = ±45°, when the Co spin is 
parallel either to the upper- or to the lower-lying oxygen 
bond; indeed, as pictorially shown in Fig. |l](a), the com- 
position of yz and zx orbitals has an asymmetric bonding 
nature in the tetrahedron, i. e. non-bonding with upper 
ligands and bonding with lower ligands or vice versa. 

Conclusions — On the grounds of the known micro- 
scopic mechanism underlying peculiar magnetoelectricity 
observed in Ba2CoGe207, we predict much stronger mag- 
netoelectric effects to appear in iron-based oxides, such 
as CaFe02, where a large polarization - magnetically con- 
trollable - is estimated. In addition to essential ingredi- 
ents, such as spin-dependent p—d hybridization and spin- 
orbit coupling, here a central role is played by the pecu- 
liar geometry, featuring "flattened" Fe04 tetrahedrons 
and non-centrosymmetric point group. The latter con- 
ditions are optimized in MgFe02, where our materials- 
design approach leads to magnified magnetoelectric ef- 
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